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Abstract: Numerous uncertainties exist in various electricity power system problems due to the size, 
complexity, geographical distribution, and influence of unforeseen events in these systems, making 
it difficult for traditional mathematics tools based on crisp set theory to have an impact on and solve 
many power system problems. As a new branch of mathematical uncertainty techniques, the 
neutrosophic expert systems approach has therefore emerged with the development of electric power 
systems and has proven successful when correctly linked. The expert typically uses ambiguous or 
neutrosophic language to describe their empirical knowledge, such as "very likely," "quite likely," "if 


wo 


x is large, then y is very likely to occur," "x should not be less than a,” etc. To design an optimal radius 


of power supply in the electrical transformer substation, this article presents a new method for 
creating primal and dual neutrosophic geometric programming problems. It also provides a 
numerical example to evaluate the approximate optimal economic power supply radius. 


Keywords: Optimal Supply Radius; Neutrosophic Geometric Programming; Neutrosophic 
Coefficients; Transformer Substations. 


NOMENCLATURES 


TS Transformer Substation 

GPP Geometric Programming Problem 

SPP Signomial Programming Problem 

NGPP_ Neutrosophic Geometric Programming Problem 

An annual extraction coefficient is the annual running cost of the entire extracted 

investment. 

a, Unconnected Part of TS's investment capacity (unit, 1Q) 

ay Investment in (11kV) Line for Every Kilometre of Unconnected Wire (unit, 1Q/km) 

by Coefficient of the Section Associated with the TS's Capacity in Investing (unit,IQ/kV A) 
Coefficient of the Wire Section-Connected Part of the (11 kV) Line Investment 
(unit, 1Q/km + mm?) 


E Coefficient of Terrain Correction 
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a Average Load Density (unit, kW /km?) 
K, Hold Ratio (the 132k V — TS takes 2.2 ~2.5) 
Pay | Mean Line Load Every Time (unit, kW) 
Uy The Middle Voltage Distribution Network's Rated Voltage (11k V) 
j Economical Current Density Wire (unit, A/mm?) 
Wire in Resistivity (unit, Q/km - mm?) 
T Annually Wasted Mean Maximum Load Hours 
Co Price of Each Wasted Watt-Hour (unit,/Q/kW.h) 
Kre Coefficient of Transformer Iron Loss(~ 0.0085kW/(kV A)?/*) 
Th Equivalent Hours of Transformer Copper Loss 
K., | Transformer Copper Loss Coefficient (~ 0.055kW/(kV A)3/*) 
Sy, | Rated Load of Transformer (unit, kV A) 
Sy Rated Capacity of Transformer (unit, kV A) 
Voltage Range in Electrically-Distributed Cities (unit, km?). 
r Mean Radius in the Supply Area, or it is the economic power supply radius 


IQ Iraqi Dinar 


1. Introduction 


The building of an optimal mathematical model for the annual cost function needs a geometric 
programming model because its variables have negative rational powers. 

It is well known that any classical geometric programming problem neglects the uncertainty part 
of the sophisticated practical environment. So, the neutrosophic model will be more flexible having 
more information. In Iraq, the services of electricity supply are very poor in almost all cities for many 
reasons, such as the cascade wars, for example, the Arab Gulf War in 1990, where the Iraqi community 
suffered from many setbacks: 

1. The economy block in the nineties of the last century, 

2. Replace the system of the government with a parliamentary system of government, 

3. Inthe twenty-one century the events of the sectarian and ethnic civil wars in Iraqi society. 

All previous reasons and more, lead to infrastructural destruction. 

This essay comes as an attempt to shed light on some important problems involving how can 
rebuild the Iraqi power electricity supply substations. The necessary point that this paper focused on 
is how to choose an optimal radius (r) of the electric power supply in the substations of Iraqi cities. 
The reader should pay attention that the monetary currency of the country is the Iraqi dinar (IQ). 

2. Classical Geometric Programming Problem (CGPP) 


To handle a class of non-linear optimization problems that are commonly seen in engineering 
design, an optimization technique known as geometric programming (GP) was developed. The GP 
technique was inspired by the work of Zener. Zener tried a novel approach in 1961 to solve a class of 
unconstrained non-linear optimization problems with polynomial terms in the objective function. To 
answer these problems, he used the well-known inequality between arithmetic and geometric means, 
which states that the arithmetic mean is greater than or equal to the geometric mean. Because the 
arithmetic-geometric mean inequality is applied, this method is referred to as the "GP technique". 
Only in cases when the objective function was unconstrained and the number of polynomial terms 
was one more than the number of variables did Zener employ this technique. Later in 1962, Duffin 
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expanded the application of this method to resolve issues when the objective function's number of 
polynomial terms is arbitrary. In 1967, Peterson enhanced the application of this technique to issues 
that also incorporate inequality restrictions expressed as posynomials with the aid of Zener and 
Duffin [1]. Additionally, Passy and Wilde (1967) extended this method to tackle issues where some 
of the posynomial terms have negative coefficients [2]. 

For instance, Duffin demonstrated that a function "duality gap" could not develop in geometric 
programming by condensing the posynomial functions to monomial form (1970) and converting 
them to linear form using a logarithmic transformation. Posynomial programming with (posy) 
monomial objective and constraint functions is equivalent to linear programming. Duffin and 
Peterson (1972) demonstrate that each of these posynomial programs GPP can be reformulated so 
that every constraint function is a (posy) binomial in that it includes at most two posynomial terms. 
An effective and extremely flexible way of solution was desired since geometric programming has 
become a popular optimization methodology. Several factors became crucial as the complexity of the 
sample geometric programs to be solved grew: Canonically, the problem's level of difficulty and 
inactive constraints [3] reported an algorithm that might take these factors into account. Later, in 1976, 
Mcnamara suggested a method for solving geometric programming problems that involved 
formulating an augmented problem with a degree of complexity zero. As a result, many algorithms 
have been proposed for solving GP, the majority of which were made before (1980); these algorithms 
are somewhat more efficient and reliable when applied to the convex problem and also avoid 
problems with derivative singularities as variables raised to fractional powers approach zero because 
logs of such variables will approach—0, There should be significant negative lower bounds on those 
variables. A significant amount of interior point (IP) algorithms for general purpose (GP) were 
developed in the 1990s. An effective method for solving posynomial geometric programming was 
developed by Rajgopal and Bricker in 2002. Condensation was a concept that was employed in the 
method, which was integrated into an algorithm for the more broad (Signomial) GP issue. The 
reformulation's constraint structure sheds light on why this algorithm is effective in avoiding every 
computing issue usually connected to dual-based algorithms. A method for addressing (positive, 
zero, or negative) variables in SPP was put out by Li and Tsai in 2005. A linear or convex relaxation 
of the original problem is computed using the majority of current global optimization methods for 
SPP. These methods might occasionally offer an impractical answer, or they might constitute the 
genuine optimum to get around these restrictions. Shen, Ma, and Chen (2008) proposed a robust 
solution algorithm for the global algorithm optimization of SPP. This algorithm adequately ensures 
the achievement of a robust optimal solution that is both feasible and close to the actual optimal 
solution and is stable under small perturbations of the constraints [4]. 

Huda E. Khalid [5] suggested an innovative GPP algorithm for discovering the ranging analysis 
by examining the impact of perturbations in the coefficients without solving the issue, as this 
proposed procedure had been caught on two coefficients at once. Additionally, the reference [6] 
investigated an original GPP employing substitution effects in a sensitivity analysis for two 
coefficients at once based on a new extended theorem and adjusted new constants. Huda E. Khalid 
[7] proposed one of the incremental procedures required to appropriately compare the results 
obtained from the incremental analysis methodology and the sensitivity analysis approach. 
Additionally, there was an effort to develop a novel computational approach that could be used to 
examine the sensitivity analysis of the geometric programming problem (GPP) and the signomial 
geometric programming problem (SPP), both of which had a difficulty level larger than zero. The 
studies [8, 9] cast a bad light on multi-objective geometric programming's degree of challenges. 
Definition 2.1: [10]: The vector x* that makes the constraint inequalities gj; (x) <1, j= 
L2iwak ,pand x >0_ into exact equalities is the optimal solution to a GPP. We refer to the constraint 
set in such issues as being tight or active. As a result, we can assess any inactive constraints by 
assessing those for which gj (x) >1 or gj; (x) <1. 
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Definition 2.2: [11] Let R%, stand for the set of real m-vectors with positive component values. Let 
there be m real positive vectors, X1,...,%m. A function with the definition f : R?, > R is called a 
aj 


monomial, and is defined as f(x) = c]]j2, xj, C> 0 where a; €R. A polynomial is a sum of 


monomials or a function of the type (x) = Miki Ce [j=1 a ,with c, >0 and a, € R. 


3. Neutrosophic Geometric programming problems ( NGPP) [4] 


In 2016, Florentin Smarandache and Huda E. Khalid developed unconstrained neutrosophic 
geometric programming, where the models were constructed as posynomials. For this subject, the 
following definitions are thought to be fundamental: 

3.1 Definition: Let h(x) be any linear or non-linear neutrosophic function, where x; € [0,1] U [0, nl] 
and x = (X1,X2,...,Xm)' an m-dimensional fuzzy neutrosophic variable vector. 

We have the inequality 

h(x) < 1 (1) 


where "<8" denotes the neutrosophic version of "<" with the linguistic interpretation being 


"less than (the original claimed), greater than (the anti-claim of the original less than), equal 
(neither the original claim nor the anti-claim)". 
The inequality (1) can be redefined as follows: 
h(x) <1 
anti (h(x)) > 1 (2) 
neut( h(x)) = 1 
3.2 Definition: Let 


N 
Xj E FN, 


The neutrosophic unconstrained posynomial geometric programming, where x = 
(X1,X2,++)Xm)* is an m-dimensional fuzzy neutrosophic variable vector, "T" represents a 
transpose symbol and h(x) = ¥}_, cx [121 x)" is a neutrosophic posynomial GP function of x , 
c, 20 a constant, yy an arbitrary real number, h(x) <¥Nz- ,\)g(x); the objective function 
h(x) can be written as a minimizing goal to consider z as an upper bound; z is an expectation 
value of the objective function h(x)," <8" denotes the neutrosophic version of "<" with the 
linguistic interpretation (see Definition 3.1), and d, > 0 denotes a flexible index of h(x). 
Note that the above program is undefined and has no solution in the case of yy <0 with some 
x;'s taking indeterminacy value, for example, 
min A(x) = 2xpexPxg? + TxPxz xz, (4) 
where x; € FN,,i = 1,2,3,4. This program is not defined at x= (.21,.3,.25,1)° , h@)= 
2(.21)~7(.3)317° + 70. 21D3(.3)75(.25) is undefined at x, =.21 with y, = —0.2. 

3.3 Definition: Let Ag be the set of all neutrosophic non-liner functions h(x) that are 
neutrosophically less than or equal to z, ie. 
Ao = {xj € FN, h(x) < ¥z}. 
The membership functions of h(x) and anti(h(x)) are: 


1 0 <h(x) <z 
h = et = st i - 
La, ( h(x)) { eds 8-2) 4 gzj(anti (@O0)-2 _ 1), gehGy <o-i4100s (5) 
0 0 <h(x) <z 


ti h = et i = a i 
Ha, (anti(h(x))) (oe (hOD)-2) _ .gg(hOo =e mid es Hie 2erd 


Eq. (6) can be changed into 
h(x) <Nz, xX = (X41, Xp, «1, Xm), Xi € FN (7) 
The above program can be redefined as follows: 
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h(x) <z 

anti(h(x)) >Z 

neut(h(x)) =Z 

X = (x1, X2, «.) Xm), Xi € FNg 

It is clear that 14,(meut(h(x))) consists of the intersection of the following functions: 


(8) 


=1 -1 t 
edo h@)-2) 21- edo Anti@))-2) (9) 
fe edo anti h@d)-2) 
-1 
eas h@-) 


z< h(x) <z—d,1n0.5 


La, (neut(h(x))) = (10) 


z—d,In0.5 < h(x) <z+d, 


3.4 Definition: Let N be a fuzzy neutrosophic set defined on [0,1] U [0, nl], € [0,1]; if there exists 
a fuzzy neutrosophic optimal point set Ag of h(x) such that 


Nw) = min{y(neut h(x))} (11) 
X = (X1,X2,«) Xm), Xi © FNg 
Nx) = edo B= RIE 2). 4 gig (@mti( Bhar CeIEx1")-2) (12) 


Then max N(x) is said to be a neutrosophic geometric programming (the unconstrained case) 
concerning N(x) of h(x). 


3.5 Definition: Let x* be an optimal solution to N(x), i.e. N(x*) = max N(x) ,x = (x1,X2,.-)Xm), Xi € 
FN,, and the fuzzy neutrosophic set N satisfying Eq. (11) is a fuzzy neutrosophic decision in Eq. (8). 


4, Designing an Optimal Radius of Power Supply in Transformers Substation Using Classical 
Geometric Programming Problems GPP [12] 


The annual-cost way had been built with consideration to the following assumptions: 

1. 132 KV-TS power supply to its consumers in a city by the direct-step-down method of 
11 KV. 

2. The load density is even over the whole electrified wire netting cover. Therefore, a static 
model is built using an annual cost as follows: 
Footy (13) 
Where Z denotes the cost of total investment, is the annual cost of investment 
operational under a certain load level, N(8—10 years) is an investment-recovery 
deadline, i.e., the total investment is returned within (8 — 10) redeemable years. 
The annual cost function in a unit capacity is denoted by 


F CoD stunts 
Ry=== N : (14) 
Where, 
Zp =a,+b,S (CQ) (15) 
Is the investment in the construction of 132 KV-TS [14]. 
Z, = L(Ma, + b2S,) (1Q) (16) 


Denotes the construction investment in the main-supply lines of the (11 KV) middle voltage 


distribution network, where M = S cos ¢ /Pyy is a circle line of (11KV) middle voltage distribution 

network, L = Er (km) is each circle-line length of it, and S; = a (mm?) denotes the wire total 
Nj 

selection in the main-supply lines of all (11 KV) middle voltage distribution net-work in (132 KV) — 

TS; S is the capacity in the TS (unit, KVA). 


MW, = H(Z, + Z;) (17) 


Ahmed K. Essa, Montifort Blessings Andrew Mitungwi, Tuweh Prince Gadama, and Ahmed A. Salama, Choosing Optimal 
Supply Radius of Transformer Substations (TSs) in Iraq’s Cities Using Geometric Programming with Neutrosophic 
Coefficients 


Neutrosophic Systems with Applications, Vol. 16, 2024 29 
An International Journal on Informatics, Decision Science, Intelligent Systems Applications 


Behaves as a direct proportion function of the unchangeable part in operations cost (large repair, 
small repair, and depreciation charge) and in the total investment of (132 KV —TS) and (11 KV) 
middle voltage distribution net line. 


= _ Epjrfo 1 (3/ -5 
Hz = APCo = 7.26 PE — 5/2 10 (18) 


Stands for the depreciation charge of (11 KV) line ina year, while by [15], the depreciation charge to 
transformers of 

(132K V—TS) is pz = (CoKpe8760 + CotpKou) X YS’ /4 NL (19) 
Where y denotes the number of transformers. When the rated capacity as Sy hours for a chosen 
transformer is 11 kV-TS, y is taken to denote an average number of transformers S/Sy in 
thel1k V —TS. 

Substitute (15)-(19) for (14), then 

Ky = (=+H)(S4+b,) +E I( +H) |— (2 ost ) #726 ee | ss (87600 Reg > 


toKe\ Pay V3 Unj Uncos?@ Vo 


3 
s°/4.NL 


CotpKeu) — (20) 


The determination of the objective function of a static (classical) model aims at making unit capacity 
annual cost minimum, ie., in Ky , with the limits to constraint S > 0, such that we have the 


following model: 
min Ky 

s.t. S>0 

Where Ky is illustrated in Eqs. (20), and (21) are called classical geometric programming models, 


(21) 


concerned with the model I. 


5. Geometric Programming with Neutrosophic Coefficients 


In this section, the objective function Kp will be reformulated from classical annual coast in unit 
capacity into minimum neutrosophic function. We expect that there are many vague or incomplete 
factors in load, investment process, and electricity prices, where they hold many neutrosophic 
phenomena. 

It is well known that the transformer substations’ capacity is non-negative (i.e. S > 0), also Ko is 
an exponential polynomial function having neutrosophic coefficients with respect to S , therefore the 
classical model can be changed into finding an answer to a neutrosophic minimum of annual cost 
under the capacity S > 0 in constraint, such that the geometric programming with neutrosophic 
coefficients can be written as a forthcoming model (22). 

5.1 Building a Neutrosophic GPP Model in an Iraqi’s Transformer Substation Depending Upon 

Neutrosophic Truth Membership Function Related to the Coefficient Cz 
Solve the following problem 


min Ko 
22 
s.t. S>0 (22) 
Where, 
Ky = 91800S-1 + 0.770795 5°5 (23) 


Having the first coefficient c, = 91800 as an ordinary real number, while the second coefficient c, = 
0.770 °° isa neutrosophic number varying in a certain interval. Suppose the truth membership 
function f4, indeterminacy membership function o,, and falsity of membership function v, for the 
coefficient c, are defined as follows: 
Let A = [0.008, 0.03] be a certain neutrosophic interval in which cz € A, 

0 


Cy < 0.008 
2 C2—0.008, 9 
Hale) =} (=) 0.008 < cy < 0.03 (24) 
1 Cy > 0.03 
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ey 0.008 < c, < 0.011 
0.011 
alco) = 42 2=004ty2 0.011 <c, < 0.03 25) 
0.008 = c, > 0.03 
1 cz < 0.008 
vs(C.) =41- a; 0.008 <c, < 0.03 (26) 
0 cz > 0.03 


It is easy to draw the above truth, indeterminacy, and falsity membership functions. 


0.5 


C2 


0.008 0.022/2 0.03 


Figure 1: The blue linear function represents the region covered by p,4(cz), while the orange function 

represents the region covered by v,(c2), and finally, the intersection region (yellow-shaded color) is 

the intersection region between p,(c,) and v4(c,) represents the region covered by o,(cz). 

On the other hand, and depending upon [16], we have 

(Sy? = 1-a > cy = 0.008 + 0.022V1—a (27) 

Here,a denotes the a — cut related to the truth membership functions f1,4(c2). 

Consequently, the program (22) is turned into 

min Ky = {91800S~* + (0.008 + 0.022V1 — @ )s°} (28) 
s.t. S>0 

Program (28) is called neutrosophic polynomial geometric programming problems (NPGPP), we 

should not forget that the degree of the dual program for eq. (28) is equal to zero, so it is immediate 

to find the solution by means of the dual program. To solve this problem more easily way, we try to 

rewrite this program (i.e. eq. 28) in its dual form as follows: 


max D = (“#0 (2.0ne40.022V1—8 oy, 
Wo 4 
S.t. Wj) +@,=1 os 


—-Wo + 0.50, = 0 

W,@, >0, a € [0,1] 
By solving the above (normality and orthogonality conditions), the values wy) = : ,O1 = “ are 
obtained. 
The well-known relationship between the primal of the GPP and its dual program is concluded by 
the below theorem: 
Theorem [13]: If (P) is canonical and there exists a x such that g,(x) <1 for = 1,2,...,p , then the 
following mathematical phrases are true: 

i. The dual program (D) has a maximizing point 6*, 
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ii. The maximum value v(6*) of problem(D) equals to the minimum value g.(x*) of (P), {ie. 
the primal and the dual problems are related through the fact that, for the optimal solution, 
min g.(x) = max v(6)}. 


iil. Each minimizing point x for (P) for the optimal solution min gg (x) = max v(6) satisfies 
6; v(6"*) for i=0 
uj(x) = 6; 


Viz=1,2,..k and 1,(6*) #0 
A, (6*) 


However, there is a good opportunity to study the various values of the dual objective function D 
depending on the values of a. The following formulas are: 
- ris the mean radius in the supply area, or it is the economic power supply radius, 


_ Sav 
r= ae (30) 


Suppose that the mean load density of Telafer township / Nineveh province/ Iraq country will raised 
to ¢ = 5195 kW/km’ by the year 2025, suppose the economic capacity of the specific transformer 
substation as S,y. 

Average item N, (which denotes the amount of 132 kV-TS that needs to be built up in Telafer district 


to meet the growing demand for electricity) given by the following formulas 

N oQKe oQKe Q (31) 
aa ra r?10Ke ~ aer42 

Keep in mind, the fact that, the area of the electrical supply on Telafer township station cover Q = 

1743.69 km?. 


- Thenumbers njina unit area in the specific transformers substation TS, 
ny = b= See (32) 
Q Ss wr; 
- Note that c, = 0.008 + 0.022V1—a@, a € [0,1], at a=0 = c, = 0.03, whileat a=1> 
c, = 0.008, as 0<a<1= 0.008 <c, < 0.03. 
The following Table 1 illustrates the maxD are the approximate optimal solutions and optimal 


values for (29), 


Table 1. The range values of D. 


Woes 
index | @ max D = (21809) ay (0-008 + 0.022/1= ai, ‘10K, 
€ [0,1] Wo W1 
K,=2.2kV; 6 
= 5915 kW/km? 
1 0 dy = 8.23119687 0.641439553 km 
2 0.1 dy = 8.02337386 0.649671351 km 
3 0.3 d3 = 7.56001993 0.66923184 km 
4 0.5 d, = 7.00574202 0.695134289 km 
5 0.7 ds = 6.29202210 0.733404644 km 
6 0.8 dg = 5.82046397 0.762466041 km 
7 0.9 d, = 5.17541713 0.808481231 km 
8 1 dg = 3.41016566 0.995585028 km 


Note that, the values of a’s were selected depending on the author’s decision. To evaluate the values 
of r;. Now we need to determine the average values S,,,, by the previous theorem. we have the 
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following relationship between the terms of the primal program and its corresponding terms of dual 


program: 
91800S-t="d; Vi=12,..,8 (33) 
2 i 
(0.008 + 0.022, /1—a)S°5 ==d,; Vi=1,2,...8 (34) 


Suppose that all evaluated S through Eq. (33) are symbolized by S,, and all computed S through 
eq. (34) symbolized by S. The following Table 2 illustrates all values of S, and S., this is for each 
i= 1,2,...8: 


Table 2. The Average Values of S; and S, Denoted by Sq, . 


_ 91800 S> 
index i Vga cP . _ 51 + Sz 
275400 2 a ote 2 
a 0.008 + 0.022,/1 — a; 
1 33458.0747 182.915482 16820.4951 
3 34324.7124 185.269297 17254.9909 
3 36428.4754 190.862453 18309.669 
4 39310.6111 198.269037 19754.4401 
5 43769.7128 209.212124 21989.4625 
6 47315.8156 217.52199 23766.6688 
7 53213.1021 230.679653 26721.8909 
8 80758.54 284.180472 40521.3603 


Substituting all amounts of S,,,, in the Table 1 to enumerate the values of 7;. Now, the last duty is to 
determine the optimal value of S,,, that gives the optimal solution for problem (28), tracking the 
concluded values of the below Table 3: 


Table 3. The approximate values of the primal program. 


index §,+5 
é wv, = : 2 min Ky = {91800(S,,;)"! + (0.008 + 0.022V1 — @ )(Sqyi)°5} 
1 16820.4951 9.34844324 
2 17254.9909 9.11264828 
3 18309.669 8.58689587 
4 19754.4401 7.95791296 
5 21989.4625 7.14789479 
23766.6688 6.61264497 
7 26721.8909 5.88038304 
8 40521.3603 3.87586527 


From Tables 2 and 3, and in terms of Eqs. (33) and (34), the approximate optimal value of Kp = 
3.8758627 which holds at the optimal economic capacity Sgyg = 40521.3603 of 132 kV- TS. 
Moreover, by the results of Table 1, we can see that the value of the approximate optimal economic 
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power supply radius is rg = 0.995585028 km, as well as the optimal value of N, = 559.96332 occurs 
at rg (ie., about (560) 132 kV —TS needs to build up in Telafer township). 


5.2 Building a Neutrosophic Model Depending Upon Either Falsity or Indeterminacy Membership 

Functions Related to the Coefficient C, 

The fact that many researchers in the field of neutrosophic optimization may ignore or may have 
been unaware of it, is that the indeterminacy membership function exactly represents the intersection 
region of the truth membership function and the falsity membership function as we sighted it by the 
representing diagram (1). Hence, for this manuscript, we took the vagueness for the second 
coefficientc,, so the following mathematical intersection is completely true: 

O4(C2) = Malco) N Valco) (35) 


For solving the program (22) with respect to its falsity membership function and depending upon the 


thoughts presented in[ J], we can suppose: 
c2—0.008\7 _ c2—0.008\2 _ = 
lo rere) =m (eo) =p (36) 
Here, B is the P- cut related to the falsity membership function v,(cz), eq. (29) is exactly given the 


same solution as Eq. (27), which means Eq. (36) is somehow considered a dual form to Eq. (27). Also, 
there is another unfathomable program that gives the solution of neutrosophic polynomial geometric 
programming (22) by using the indeterminacy membership function o,(c,): 


1 c2—0.008\7 _ 
a 0.022 ) =r : (37) 
s 2 F —0.008 1 
Where y € [0,1] is the y- cut corresponding to o,(.), Eq. (30) implies that (2 ) SS SG 
0.011 = 0.0009,/0.5-—y => 
Cp = 0.011 + 0.0009,/0.5 — y. 
Consequently, the program (15) can be reformulated as: 
min Ky = {91800S~! + (0.011 + 0.0009,/0.5 — y )S°5} Ae 


s.t. S>0 
Program (38) is defined as a neutrosophic posynomial geometric programming problem, and the 
main difference between program (28) and program (38): is that program (28) is NPGP related to the 
truth-membership function of the coefficientc,, while the program (38) is NPGP related to the 
indeterminacy membership function of the same coefficientc,. Also, the same previous analysis tables 
(1, 2, 3) that were for a,maxD, r; can be re-written and re-analysis for y, max D ,7;. 


6. Conclusions 


This manuscript discussed an innovative trying to analyze the annual cost of investment in the 
power supply systems with respect to the validity of some uncertainty by assuming the neutrosophic 
coefficient cz, to be able to calculate the mean of the economic capacity and economic supply radius 
of 132 kV- transformer substation in Telafer township, all these issues held by building mathematical 
non-linear programming named Neutrosophic Geometric Programming Problems (NGPP) this was 
for the first time, also there more analyzing techniques can be made for the same geometric 
programming problems with the neutrosophic point of view, either by a truth membership function, 
or by a falsity membership function, or by an indeterminacy membership function, to their 
corresponding a, B,y- cuts. 
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